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ABSTRACT 

The integrated line width derived from CO spectroscopy provides a powerful tool to study the internal kine- 
matics of extragalactic objects, including quasars at high redshift, provided that the observed line width can be 
properly translated to more conventionally used kinematical parameters of galaxies. We show, through construc- 
tion of a /Tj-band CO Tully-Fisher relation for nearby galaxies spanning a wide range in infrared luminosity, 
that the CO line width measured at 20% of the peak intensity, when corrected for inclination and other effects, 
successfully recovers the maximum rotation velocity of the disk. The line width at 50% of the peak intensity 
performs much more poorly, in large part because CO lines have a wide range of profiles, which are shown to 
vary systematically with infrared luminosity. We present a practical prescription for converting observed CO line 
widths into the stellar velocity dispersion of the bulge (cr*), and then apply it to a sample of low-redshift iz <, 
0.2) and high-redshift (1-4 < z 6.4) quasars to study their host galaxies. Nearby quasars roughly fall on the 
correlation between black hole mass and bulge stellar velocity dispersion established for inactive galaxies, but the 
host galaxies of the high-z quasars systematically deviate from the local M, -cr* relation. At a given cr*, high-z 
quasars have black hole masses larger by a factor of ~4 relative to local galaxies, suggesting that early in the 
life-cycle of galaxies the development of the bulge lags behind the growth of the central black hole. An alternative 
explanation for these observations, which currently cannot be ruled out rigorously, is that high-redshift quasars 
are preferentially viewed at face-on orientations. 

Subject headings: galaxies: bulges — galaxies: ISM — galaxies: kinematics and dynamics — galaxies: nuclei 
— (galaxies:) quasars: general — galaxies: Seyfert 
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L INTRODUCTION 

The discovery of strong empirical relations between black 
hole (BH) mass and the properties of nearby galaxies, both in- 
active (Kormendy & Richstone 1995; Magorrian et al. 1998; 
Gebhardt et al. 2000a; Ferrarese & Merritt 2000) and active 
(Gebhardt et al. 2000b; Ferrarese et al. 2001; Barth et al. 
2005; Greene & Ho 2006), has triggered an upsurge of inter- 
est in the coevolution of BHs and galaxies (e.g.. Ho 2004). An 
important unanswered question is when the local scaling re- 
lations were established. While there is no shortage of opin- 
ions from theoretical studies (e.g., Croton 2006; Fontanot et 
al. 2006; Hopkins et al. 2007), to date solid empirical evi- 
dence has been relatively meager, and somewhat conflicting. 
From stellar velocity dispersion measurements of moderate- 
luminosity active galactic nuclei (AGNs), Treu et al. (2004, 
2007) and Woo et al. (2006) find that the M, - cr* relation at 
z = 0.36 already begins to show departure from the local rela- 
tion. On the other hand, studies of larger and more luminous 
AGN samples, but using less direct surrogates of a^, based on 
the widths of optical nebular emission lines, find little or no 
evolution of the M, - cr* relation out to z « 1.2 (Salviander 
et al. 2007) and z « 3 (Shields et al. 2003). At even higher 
redshifts, Walter et al. (2004) resolved the molecular gas in 
the z = 6.42 quasai- SDSS J114816.64H-525150.3 into an appar- 
ent nuclear disk, from which they were able to derive a dy- 
namical mass that, taken at face value, indicates that the host 
galaxy has an embryonic bulge compared to the estimated mass 
of the BH. WeiB et al. (2007) arrived at a very similar con- 
clusion from their detailed study of the z = 3.91 lensed quasar 
APM 08279H-5255. While the interpretation of these results is 
not unambiguous, it qualitatively agrees with the photometric 



and lensing studies of high-z quasars by Peng et al. (2006a, 
2006b), which show that at 1 z 4.5 the ratio of BH mass to 
host galaxy mass is a factor of several higher than the value in 
the local Universe. 

The masses of high-z quasar host galaxies are challenging to 
estimate accurately, and likely will remain so for some time. 
To make progress one should explore multiple avenues. One 
promising technique is to exploit the kinematic information en- 
coded in molecular line observations. While studying molec- 
ular gas in high-z sources is still far from routine, an increas- 
ing number of detections have been reported (Solomon & Van- 
den Bout 2005, and references therein), and the upcoming gen- 
eration of new millimeter-wave facilities hold great promise 
for rapid progress in this area. Shields et al. (2006) recently 
attempted to constrain the host galaxy masses of nine high-z 
quasars with available measurements of the CO line width and 
BH mass. In order to place the high-z objects on the M, - a^, 
relation, they calibrated the width of the CO line profile at 50% 
of the peak intensity (hereafter W50) against a small sample of 
low-z quasars for which they could obtain indirect estimates of 
the bulge stellar velocity dispersion, concluding in the end that 
(7* could be estimated simply from W50/2.35. Since the cur- 
rently available high-z quasars have conspicuously narrow CO 
lines (Greve et al. 2005; Carilli & Wang 2006), Shields et al. 
(2006) suggest that these objects deviate strongly from the local 
M, - (T, relation in having a much larger BH mass for a given 
(7* . The implication is that the growth of the BH predates, and 
is largely decoupled from, the formation of the galaxy bulge. 

Although the interpretation of Shields et al. (2006) res- 
onates well with the findings of Peng et al. (2006a, 2006b), 
it is not without complications. With a few exceptions (e.g., 
SDSS J114816.64H-525150.3; Walter et al. 2004), most of 
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the other quasars are not spatially resolved, and thus we have 
no knowledge of the spatial distribution of the molecular gas. 
Apart from signifying a shallow gravitational potential, a nar- 
row CO line can arise from at least two other effects. First, the 
CO may be sufficiently concentrated toward the center of the 
host that it only samples the innermost part of the galaxy's ro- 
tation curve, thereby underestimating the maximum rotation ve- 
locity. This is well-known from observations (e.g., Downes & 
Solomon 1998) as well as numerical simulations (e.g., Barnes 
& Hernquist 1996) of merging galaxies, wherein the bulk of 
the molecular gas, as a result of tidal torques, sinks to the cen- 
tral few hundred parsecs and forms a compact nuclear disk or 
ring. As the hosts of high-z quasars likely involve mergers, this 
should be a source of concern. Second, as argued by Carilli 
& Wang (2006), the distribution of inclination angles for the 
CO-emitting disk of high-z quasars, by virtue of their luminos- 
ity bias and selection method, may preferentially favor more 
face-on orientations. If so, then the inclination correction to the 
CO line widths can be arbitrarily large. For the high-z sample 
to conform to the local M, - cr, relation, the average inclina- 
tion angle needs to be (/) « 13-15° (CarilU & Wang 2006; Wu 
2007). Lastly, we note that the interpretation of these and future 
observations depends quite critically on how one actually trans- 
lates the CO line width into a more familiar dynamical variable. 
Shields et al. (2006) chose to calibrate W50 with three indi- 
rect estimators of cr* (virial BH masses and [O III] A5007 line 
widths for seven objects, and host galaxy luminosities for two 
objects). By contrast, Wu (2007) examined directly the correla- 
tion between W50 and cr* for a significantly larger sample of 33 
Seyfert galaxies, concluding that the Shields et al. calibration 
is significantly nonlinear 

This paper presents a new method of estimating cr* based on 
more complete information on the CO line profile. We begin 
by asking three basic questions that we believe have not been 
adequately addressed. (1) How well does the CO line width 
trace the kinematics of galaxies and precisely what part of the 
galaxy does it probe? (2) What part of the CO profile, which, 
unlike that of H I is often not sharp and double-horned, should 
one measure? And (3) how do we relate the CO line width to 
the kinematics of the bulge? We investigate the first two issues 
by reexamining the CO Tully-Fisher relation, with particular 
emphasis on infrared-luminous galaxies, perhaps the best local 
analogs of what high-z quasar host galaxies might look like. We 
find that galaxies over a wide range of infrared (IR) luminosities 
fall on the CO Tully-Fisher relation provided that one uses W20, 
the width of the CO line at 20% of the peak intensity, rather 
than Wso. Moreover, we show that the shape of the line profile, 
as parameterized by the ratio W20/W50, varies significantly and 
systematically with IR luminosity. This thus provides a sim- 
ple prescription for converting W50 to the preferred parameter 
1^20 if the IR luminosity is known. Having shown that the CO 
line width can effectively probe the maximum rotation veloc- 
ity (Vm) of the galaxy, we then use the new calibration of the 
Vm-c* relation by Ho (2007) and discuss the implications for 
high-z quasar host galaxies. 
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INFRARED-LUMINOUS GALAXIES 

The CO Tully-Fisher relation is analogous to the more famil- 
iar relation commonly used with H I (TuUy & Fisher 1977). 
Originally proposed by Dickey & Kases (1992) and Sofue 
(1992), it has been applied to normal galaxies both nearby 
(Schoniger & Sofue 1994, 1997; Lavezzi & Dickey 1998) and 
at moderate distances up to z ~ 0.1 (Sofue et al. 1996; Tutui 
et al. 2001). As discussed by these authors, CO offers a num- 
ber of potential advantages compared to H I, chief among them 
being the relative insensitivity of the CO distribution to environ- 
ment perturbations and the technical feasibility of observing to 
greater distances. 

One of the practical difficulties of using the CO line stems 
from the wide variety of line profiles actually observed. Unlike 
H I, which generally exhibits a sharp-edged, double-horned pro- 
file characteristic of an optically thin rotating disk, CO line pro- 
files are much more irregular, with shapes ranging from classic 
double peaks, to more rounded single peaks, to near Gaussians 
(e.g., Sanders et al. 1991; Dickey & Kases 1992). As discussed 
in Lavezzi & Dickey (1997), many factors contribute to this 
diversity in CO line profiles, two physically interesting ones 
being the radial distribution of the gas and its opacity. In spi- 
ral galaxies the molecular gas is generally more centrally con- 
centrated than the atomic hydrogen, and the '^CO(l-O) transi- 
tion tends to be optically thick. Nevertheless, in practice the 
full-width near zero intensity (or approximately W20) of the CO 
line quite faithfully follows that of H I (Dickey & Kases 1992; 
Schoniger & Sofue 1994). The molecular gas distribution, al- 
though more concentrated than the H I, evidently extends suffi- 
ciently past the peak of the rotation curve to imprint the max- 
imum rotation velocity on the CO line profile. To achieve the 
best possible match between CO and H I velocities, one should 
apply corrections for interstellar turbulence (Lavezzi & Dickey 
1997), as well as accounting for the relative differences be- 
tween the radial distribution of the two tracers and luminosity- 
dependent differences in the gradient of the rotation curve (Tu- 
tui & Sofue 1999). 

Do quasar host galaxies also obey the CO Tully-Fisher re- 
lation? Insofar as galaxy interactions or mergers are often im- 
plicated as the triggering mechanism of quasars (e.g., Hopkins 
et al. 2006), which are accompanied or preceded by starburst 
activity (e.g.. Ho 2005), this issue can be effectively addressed 
by looking into the properties of nearby IR-luminous galaxies 
(Sanders & Mirabel 1996). We have assembled a large sam- 
ple of galaxies spanning a wide range in IR luminositie^ that 
contain measurements of both W20 and W50 for the '^CO(l- 
0) line integrated over at least a central diameter of ~10 kpc, 
but preferably covering the entire galaxy. This last criterion is 
very important because most nearby galaxies are larger than the 
beam sizes of millimeter-wave telescopes, whereas we wish to 
simulate, as closely as possible, the global measurements typ- 
ically achieved with single-dish H I telescopes in order to cap- 
ture the full extent of the rotation curve. The rotation curves of 
spiral galaxies usually reach their maximum velocity at r « 2-7 
kpc, beyond which they remain flat (Sofue & Rubin 2001). To 
achieve this, most of the objects are relatively distant (60% are 
beyond 70 Mpc); more nearby systems have been mapped, and 
the published line widths represent the synthesized profile over 
the mapped positions. The final sample of 151 galaxies, assem- 



'Throughout this paper Ljr refers to the total IR (8-1000 ^m) luminosity, derived using 12, 25, 60, and 100 im\ flux densities from the Infrared Astro- 
nomical Satellite (IRAS) following the prescription of Sanders & Mirabel (1996). The IRAS data are taken from the NASA/IPAC Extragalactic Database (NED; 
http : / / nedwww . ipac . caltech . eduj. 
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Note. — Col. (1) Galaxy name. Col. (2) Luminosity distance, derived from the mean heliocentric velocity listed in 
HyperLeda (http : //leda . univ-lyonl . f r) assuming a Local Group infall velocity of 208 km s~^ toward the Virgo cluster 
and a Hubble constant of Hq = 70 km s~^ Mpc~^. Col. (3) Inclination angle. Entries followed by a colon are especially 
uncertain because the galaxy morphology is severely disturbed. The inclination angles for the objects from Lavezzi & 
Dickey 1998 are not listed, but they were applied by these authors to the line widths in Cols. (4) and (5). Col. (4) CO line 
width at 20% of the peak intensity, corrected for inclination angle and interstellar turbulence following the prescription 
of Lavazzi & Dickey 1997. Note that these values of W^lo have not yet been corrected for the velocity-dependent effect 
discussed by Tutui & Sofue (1999; see Eq. 1). Col. (5) Same as for Col. (4), but for the CO line width at 50% of the peak 
intensity. Col. (6) Total absolute -ft's-band magnitude from 2MASS; these values have not been corrected for Galactic or 
internal extinction. Col. (7) Total IR (8-1000 jim) luminosity. Col. (8) References for the CO data; (1) Sanders et al. 
1991; (2) Schoniger & Sofue 1994; (3) Lavazzi & Dickey 1998. Col. (9) Notes: (a) Extreme low-velocity outher in the 
TuUy-Fisher relation; (b) IRAS 25 /im upper limit; (c) IRAS 12 fim upper limit; (d) K-hand data from Kim et al. 2002; 
(e) Inclination angle from McLeod & Rieke 1995; (f) Inclination angle from Scharwachter et al. 2003; (g) i^-band data 
from Bushouse & Stanford 1992; (li) A'-l)aiid data from McLeod & Rieke 1995; (i) A'-band data from Busliouse & Werner 
1990. Table 1 is presented in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. 



bled in Table 1, covers Lir w 10^ to 4 x 10^^ L©, ranging from 
nearby normal galaxies to those of the "luminous IR" and "ul- 
traluminous IR" variety. The data were taken from three compi- 
lations: (1) luminous //JA^ galaxies from Sanders etal. (1991); 
(2) nearby galaxies from Schoniger & Sofue (1994), which in 
addition contains the two powerful AGNs Mrk 23 1 and I Zw 
1; and (3) the distant cluster members and ultraluminous IR 
galaxies from Lavezzi & Dickey (1998). 

Inclination angles were estimated from axial ratios obtained 
mostly from de Vaucouleurs et al. (1991), and otherwise from 
Paturel et al. (2000), as described in Ho (2007), except that the 
intrinsic thickness of the disk was fixed at a value of 0.2 instead 
of allowing it to vary with Hubble type. The morphologies of 
many of the IR-luminous galaxies from the sample of Sanders 
et al. (1991) are particularly difficult to define, as a number 
of them are merger remnants or otherwise appear tidally dis- 
turbed. In Table 1 we have flagged the most suspicious values 
with a colon. To avoid excessively large and uncertain inch- 
nation corrections, we excluded objects with inclination angles 
less than 30°, except in the case of the Sanders et al. (1991) 
sample on account of the inherent difficulty of obtaining accu- 
rate inclinations for many of these sources. Following Lavezzi 
& Dickey (1997), we apply a correction for turbulence broad- 
ening to the CO line widths, adopting their formalism and their 
recommended value of 38 km s"' and 19 km s"' for the cor- 
rection to W20 and W50, respectively. [We did not apply correc- 
tions for inclination or turbulence broadening to the Lavezzi & 
Dickey (1998) sample because these authors had already done 
so.] Lastly, we also apply the velocity-dependent correction 
recommended by Tutui & Sofue (1999) to make W20 fully com- 
patible with the H 1-based velocities. To account for the differ- 



ent radial distribution of CO versus H I and for variations in the 
gradient of the rotation curve with luminosity (and hence rota- 
tion velocity), these authors suggest the following cahbration: 

1^20 = 0-76 1^20 + 83. 8. (1) 

We confirm that applying this additional correction does result 
in marked improvement in the final CO TuUy-Fisher relation. 

Figure 1 shows the /Tj-band Tully-Fisher relation for the sam- 
ple, plotting separately the maximum rotation velocity as es- 
timated from the corrected CO line widths measured at 50% 
(W^q) and 20% (W20) of the Une peak, and dividing the sam- 
ple into normal and IR-luminous galaxies. With just a few ex- 
ceptions, the Ks-hand (2.16 fim) magnitudes come from homo- 
geneous photometry from the Extended Source Catalog of the 
Two-Micron All-Sky Survey (2MASS; Skrutskie et al. 2006). 
Overplotted on the figure is the H I TuUy-Fisher relation derived 
by Verheijen (2001) for galaxies in the Ursa Major cluster in the 
K' band, which is quite similar to the 2MASS band (Bessel 
2005). Looking first at the normal galaxies (top two panels), it 
is clear that fares better than W^q in matching Verheijen's 
relation. Not only does the zeropoint agree better, the scatter 
is also reduced, from a root-mean-square (rms) of 1 .0 mag to 
0.77 mag. A large part of the improvement stems from the fact 
that the Tutui-Sofue correction has been applied to W20 but not 
to W50 (because it is not currently known for the latter); without 
this correction, the scatter reduces only to an rms of 0.93 mag. 
The observed scatter for our sample is larger than Verheijen's, 
but this is hardly surprising considering the larger distance er- 
rors and greater heterogeneity of our sample. The bulk of the 
objects lie within the boundaries that enclose twice the rms 
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Fig. 1 . — The CO TuUy-Fisher relation in the Ks band for (top) normal galaxies and (bottom) IR-luminous and IR-ultraluminous galaxies. The maximum rotation 
velocity is estimated using {left) i Wjq and {right) i . The values of W^q and have been corrected for inclination and turbulence, and has been corrected 
in addition for the Tutui-Sofue (1999) effect (see Eq. 1). In the top panels the circles are the nearby galaxies from Schoniger & Sofue (1994), and the squares are 
the cluster galaxies from Lavezzi & Dickey (1998); only objects with i > 30° are included. In the bottom panels, the triangles come from the sample of Sanders 
et al. (1991), where the solid symbols represent objects without obvious optical morphological peculiarities, open symbols represent those with clearly disturbed 
morphologies or tidal features, and skeletal symbols denote those with inclination angles formally less than 30°. The open stars mark the 21 ultraluminous IR 
galaxies from the study of Solomon et al. (1997) for which Lavezzi & Dickey (1998) obtained follow-up optical imaging to determine inclination angles. Mrk 231 
and I Zw 1 from Schoniger & Sofue (1994) are plotted as an asterisk and a plus symbol, respectively. The solid line shows the K'-haai TuUy-Fisher relation from 
Verheijen (2001) based on H I rotation curves of spiral galaxies in the Ursa Major cluster; the dashed lines lines mark the region that has twice the rms scatter of the 
Ursa Major objects. 



scatter of Verheijen's "H I" sample (0.59 mag). There appears 
to be a small zeropoint offset of ~ 0.3 mag (or ~ 9% increase in 
velocity), but given the limited size and character of our sample, 
we regard this finding as tentative. Overall, our analysis reaf- 
firms the results of previous attempts to compare the CO and 
H I TuUy-Fisher relations at optical wavelengths, now extended 



for the first time to the near-lR. 

That 1^20 superior than W^q in recovering the maximum 
rotation velocity is most dramatically illustrated for the sam- 
ple of IR-luminous galaxies shown on the bottom two panels of 
Figure 1 . Employing W^q places the bulk of the objects signif- 
icantly offset (by ~ 2 mag), with enormous scatter (rms =1.9 
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mag), to the low-velocity side of the reference Tully-Fisher lo- 
cus. The use of Wjg, on the other hand, shifts most of the objects 
back to the fiducial ridgeline and markedly tightens the distribu- 
tion (rms = 1.5 mag). Again, the Tutui-Sofue correction to W20 
has a significant impact; had it not been applied, the rms scatter 
is reduced only to 1.7 mag. To be sure, the scatter is larger than 
in the case of the normal galaxies, but this is entirely to be ex- 
pected given how poorly we know the inclination angles, not to 
mention of the strong tidal perturbations witnessed in the stars, 
which presumably must also be imprinted to some degree on 
the gas. There are no clear-cut differences when the sample is 
separated either by the degree of morphological peculiarity as 
judged by our qualitative visual examination of optical images 
in NED (open versus solid triangles) or by the more detailed 
"interaction types" given in Sanders et al. (1991), nor do the 
highest luminosity sources (open stars) stand out in any con- 
spicuous manner. Despite the appreciable star formation rates 
in many of these systems, most of which qualify as hefty star- 
bursts, we see no luminosity offset relative to the fiducial ridge- 
line. The modest amount of brightening predicted in the K band 
for starbursts (~ 0.5 mag; Mouhcine & Langon 2003) may be 
imperceptible in the midst of the large scatter in the figure. 

As an aside, we note that our sample of normal galaxies 
(Fig. lb) contains 3-4 outliers with exceptionally low rota- 
tion velocities, reminiscent of the population of kinematically 
anomalous galaxies identified by Ho (2007) on the basis of H I 
measurements. Ho posits that these objects possess substantial 
amounts of dynamically unrelaxed atomic hydrogen captured 
either from minor mergers or possibly accreted directly from 
primordial clouds. However, in view of the substantial quan- 
tity of molecular gas involved in the current sample, a primor- 
dial origin seems far-fetched. Instead, we may appeal to the 
other arguments presented by Ho (2007), namely that the gas 
and stars might be strongly misaUgned, or else the gas has been 
severely perturbed by tidal interactions. Since all of the outliers 
belong to the cluster sample of Lavezzi & Dickey (1998), this is 
perhaps not implausible. The IR-luminous subsample (Fig. Id) 
also contains a number of low-velocity outliers, but these are 
balanced by an almost equal number of high- velocity outliers, 
a situation that can be explained by the gas and stars not be- 
ing coplanar (Ho 2007). Given the great difficulty in estimating 
inclination angles for this subsample, however, we should not 
attach too much weight to the observed scatter. 

In the context of this study, the critical point of the above 
analysis is that the CO line width, as measured by Wjq, pro- 
vides a robust measurement of the maximum rotation velocity 
of the large-scale disk of the galaxy, on nearly equal footing 
with velocities more conventionally derived from integrated H I 
Une widths or extended H I or Ha rotation curves. This holds 
not only for normal galaxies, but also for IR-luminous systems, 
whose evolutionary status may be hoped to bear some resem- 
blance to that of quasar host galaxies. 

3. ESTIMATING a» FROM CO LINE WIDTHS 

We have established that the width of the CO line measured 
at 20% of the peak intensity, W2Q, accurately recovers the max- 
imum rotation velocity of the disk, On the other hand, the 
kinematic quantity most closely linked to BH mass is the ve- 
locity dispersion of the bulge, cr*, not Vm- A number of stud- 
ies (e.g., Ferrarese 2002; Pizzella et al. 2005) have suggested 
that galaxies over a wide range of morphological types obey 
a tight, nearly linear correlation between ti„ and ct*, and by 



virtue of the M, - cr* relation, that BH mass may be equally, 
or perhaps even more fundamentally, coupled to than to cr, . 
The recent analysis of Ho (2007; see also Courteau et al. 2007), 
however, shows that in fact the v^-o'* relation exhibits substan- 
tial intrinsic scatter and that its zeropoint varies systematically 
with galaxy morphology, bulge-to-disk ratio, and light concen- 
tration. Despite these complications, correlated with a*, 
both as observed, not only for inactive but also active galax- 
ies (Ho et al. 2007), and as anticipated from basic dynamical 
considerations (see discussion in Ho 2007), and so we can take 
advantage of this fact to translate the CO line width into tr*. 
Ho (2007) presented fits to the Vm-ao relation using the largest 
and most comprehensive calibration sample to date. Since the 
zeropoint depends on morphological type, and the scatter is 
somewhat reduced for later-type galaxies, we use here the fit 
for the 550 "kinematically normal" spiral galaxies in Ho's sam- 
ple (Eq. 5 in Ho 2007), to yield 

logo-* = (1.26±0.046)logv„-(0.78±0.11), (2) 

where in the present context v„, w jW-^q- 

Despite the advantages of over W^q, two practical diffi- 
culties often preclude using the former over the latter. The line 
width near the base of the profile is maximally sensitive to the 
noise in the spectrum and to errors in baseline determination, 
and hence W20 is inherently more uncertain than W50. Moreover, 
many studies, presumably for this very reason, do not even pub- 
lish W20, but W50 is often given. Now, 1^20 can be estimated from 
1^50 in two simple cases: for a sharp-edged, double-horned pro- 
file W20 /W50 « 1 , and for a Gaussian W20 / W50 « 1 .5 . Many CO 
profiles, however, do not conform to these two simple shapes, 
and indeed real line profiles can even have wings more extended 
than a Gaussian, and so it would be highly desirable to devise 
an empirical guideline to estimate W20 from ^^50. Note that the 
differences between these two examples have dramatic conse- 
quences for the inferred BH masses. A 50% increase in the 
CO Une width, for example, when folded through Equation 2 
and the M. - tr* relation of Tremaine et al. (2002; M. oc ct^ *'^), 
results in a factor of ~ 8 increase in the predicted BH mass. 

Recall from our discussion of Figure 1 that the impact of 
substituting W^q with had a much more dramatic effect in 
recovering the H I Tully-Fisher relation in the case of the IR- 
luminous galaxies than for normal galaxies. The implication is 
that the line shape must vary systematically with IR luminos- 
ity. We illustrate this effect in Figure 2, where we plot WjqAVjq 
as a function of the total IR luminosity. At Ljr ^ lO'" Lq, 
Wjq/Wjq ~ 1, but W20AV5Q progressively increases toward higher 
luminosities, such that by Lir « 10'^ Lq it attains a median 
value close to 1.5, the expected value for a Gaussian adopted 
by Shields et al. (2006). Despite the large scatter and the fact 
that ^20^^50 is bounded on one side, a generalized Spearman 
rank-order correlation test formally yields a correlation coeffi- 
cient of p= 0.48, with a probability of < 10"** for rejecting the 
null hypothesis of no correlation. A linear regression fit with 
logLiR as the independent variable gives 

W^2oMo = (0.21 ±0.041)logLiR-(1.05 ±0.44). (3) 

Note that at any given luminosity there are objects spanning the 
full range of Wjq/Wjq, from ~ 1 up to a maximum value that in- 
creases roughly monotonically with luminosity. It is the upper 
envelope of the Wjo/Wjo distribution, and consequently the 
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9 10 11 12 

log (^IR / Lq) 

Fig. 2. — Variation of the line sfiape, as parameterized by tlie ratio WjqAVjq, 
as a function of total IR luminosity. Tlie symbols have the same meaning as in 
Figure 1. The solid line shows the linear regression fit between W^q/W^q and 
logLiR (Eq. 3). 

dispersion about the mean, that increases with luminosity; at 
the highest luminosities ^20^50 can be as large as 2-2.5. Since 
physically Wjq should be larger than W^q, even if measurement 
errors may result in Wjq/Wjq ^ 1, Equation 3 strictly speaking 
should only be applied when ^2^0/^50 > 1- Also, it is worth 
stressing that the correction recommended in Equation 3 is sta- 
tistical in nature; given the large scatter about the mean relation, 
especially at high luminosities, applying the correction to any 
individual source is subject to a large uncertainty. 

Interestingly, Lavezzi & Dickey (1997) searched for possi- 
ble correlations of CO line shape with galaxy properties but 
found none that was especially compelling. It could be that 
the effect described here was masked by their smaller sample 
and more limited dynamic range. Multiple factors influence 
the observed CO line profile (Dickey & Kases 1992; Lavezzi 
& Dickey 1997), but chief among them are the optical depth 
and the radial distribution of the gas. The optical depth of the 
'^CO(l-O) line in starbursts, however, if anything, appears to 
be lower than in normal galaxies (e.g., Aalto et al. 1995; Glenn 
& Hunter 2001), in which case Figure 2 actually exhibits the 
opposite trend predicted if opacity were the principal cause for 
the Une shape variation. On the other hand, it is entirely natural 
to expect the molecular gas to be more centrally concentrated in 
starburst galaxies. The most powerful starbursts seem to arise 
from central stockpiles of molecular gas that have been dissi- 
pated to the nuclear region, often with the aid of tidal interac- 
tions (e.g., Barnes & Hemquist 1996; Sanders & Mirabel 1996; 
Downes & Solomon 1998). 

In summary, we recommend the following recipe for con- 
verting an observed CO line width into a*. This recipe is em- 
pirically motivated by the success with which we were able to 
recover the rotational velocity of IR-luminous galaxies, likely 
good analogs of quasars, as described in Section 2. If the line 
width at 20% of the peak can be reliably measured, then W20 
should be used, after correcting it for turbulence as described 
in Lavezzi & Dickey (1997), inchnation, and the Tutui-Sofue 
(1999) effect (Eq. 1). If only the width at 50% of the peak is 
known or can be measured with confidence, then after correct- 
ing for turbulence and inclination. Equation 3 should be used in 
concert with an estimate of Ljr to convert W^q to Wjq, followed 
by an adjustment for the Tutui-Sofue (1999) effect. If the incli- 



nation angle of the source is unknown, as is almost always the 
case in high-z objects, but can be assumed to be random, then 
(i) =45°. With 11^20 hand, ct* follows from the v^-a^ relation 
(Eq. 2). 

4. THE HOST GALAXIES OF HIGH-Z QUASARS 

We apply the recipe outlined above to reassess the M, - cr* 
relation of high-z quasars. For comparison, we also evaluate all 
the low-z quasars with available CO line widths. The data, com- 
piled in Table 2, are similar to those presented in Shields et al. 
(2006), but with several important modifications. First, we ex- 
panded the sample of low-z Palomar-Green (PG) quasars from 
12 to 17. The five additional sources have uniformly measured 
broad H/? Une widths from Boroson & Green (1992) and optical 
spectrophotometry from Neugebauer et al. (1987), and so BH 
masses can be derived for them using the conventional virial 
method. The line widths from the original sample of Evans et 
al. (2001) have now been updated by Evans et al. (2006). For 
ease of comparison, we adopt the particular virial relation used 
by Shields et al., as well as their choice of cosmological param- 
eters to derive distances. One key difference lies in the treat- 
ment of the CO line widths. The line widths for the nine ob- 
jects published in Scoville et al. (2003) pertain to the full-width 
near zero intensity (~ W20). not the full-width at half maximum 
(W50). Shields et al. obtained the Wso values of these sources 
by assuming that W50 = |W2o. Since the formalism developed 
in this paper explicitly accounts for line shape variations. Ta- 
ble 2 carefully documents the original data. We also found it 
necessary to adjust the line width data for a few of the high-z 
quasars. 

Figure 3 summarizes the CO-derived M, - a* relation for 
quasars, where we have appUed a correction for an average in- 
clination angle of 45°. The low-z PG quasars roughly center 
around the local M, - cr, relation of inactive galaxies (Tremaine 
et al. 2002), but the scatter is very large. Taken at face value, 
the PG quasars actually he slightly below the local relation, 
reminiscent of a similar tendency reported in other broad-Une 
AGNs (Greene & Ho 2006). We refrain from quantifying the 
exact magnitude of the offset, as it depends on the choice of 
the virial relation used to obtain the BH mass, and because 
the current sample is both small and compromised by the lack 
of inchnation angle information. Two of the objects have ex- 
ceptionally narrow CO lines (PG 0838H-770, W50 = 60 km s"'; 
PG 1415H-451, W50 = 90km s~'), but they should not be viewed 
as too alarming, since similar cases are seen among inactive 
normal and IR-luminous galaxies. Both of these quasars pos- 
sess luminous hosts (Mk « —25.8 and —24.7 mag, respectively, 
assuming H-K k, 0.3 mag; McLeod & Rieke 1994), which 
would place them among the kinematically anomalous mem- 
bers in the CO TuUy-Fisher relation (Fig. lb). 

More interesting is the situation for the high-z quasars. Con- 
sistent with the analysis of Shields et al. (2006), high-z quasars 
seem to show a clear statistical offset above the local M, - cr* 
relation. Our treatment of the CO hue widths, however, has 
reduced considerably the magnitude of the offset, from an av- 
erage value of AlogMfiH =1.8 dex as found by Shields et al. 
to AlogMfiH ~ 0.6 dex. This brings these results into better 
agreement with the findings of Peng et al. (2006b), whose pho- 
tometric analysis of high-z quasar host galaxies suggests that at 
1.7 < z < 4.5 the ratio of BH mass to bulge stellar mass is a 
factor of 4-6 higher than the local value. Of course, without 
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TABLE 2: The Quasar Sample 



Quasar 


z 


log Mbh 


^^20 








lOgilR 


Reference 


Notes 








(km s-i) 


(km s^"*") 


(km s-i) 


(km s^"*^) 


(i0) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


PG 0050+124 


0.0611 


6.95 


405 


400 


596 


218 


11.93 


1 


a 


PG 0052+251 


0.171 


8.61 


500 




580 


210 


11.72 


3 


b 


PG 0157+001 


0.163 


8.02 




270 


517 


182 


12.65 


4 




PG 0804+761 


0.100 


8.02 


881 




990 


412 


11.72 


2 




PG 0838+770 


0.131 


7.71 




60 


164 


43 


11.60 


4 




PG 1119+120 


0.0502 


7.22 




220 


364 


117 


11.16 


4 




PG 1229+204 


0.0630 


7.87 


266 




329 


103 


11.26 


2 




PG 1351+640 


0.088 


8.67 


414 




488 


169 


11.87 


2 




PG 1404+226 


0.0980 


6.71 


468 




546 


195 


11.12 


2 




PG 1415+451 


0.114 


7.68 




90 


203 


56 


11.50 


4 




PG 1426+015 


0.0865 


8.63 


633 




723 


278 


11.51 


2 




PG 1440+356 


0.0791 


7.18 


600 




688 


261 


11.57 


2 




PG 1444+407 


0.267 


8.14 


300 




366 


118 


12.34 


3 


b 


PG 1613+658 


0.129 


8.94 




400 


688 


260 


12.02 


4 




PG 2130+099 


0.0629 


7.61 


487 




566 


204 


11.45 


2 




PG 2214+139 


0.0658 


8.25 


488 




567 


205 


11.08 


2 




Q 0957+561 


1.412 


8.86 


660 




752 


292 


12.81 


5 


c 


Cloverleaf 


2.558 


8.2 




416 


777 


304 


12.73 


6 


d 


RXJ J0911+0551 


2.796 


8.85 




350 


634 


235 


12.36 


7 




SMM J04135+1027 


2.846 


8.82 




340 


691 


262 


13.38 


7 




APM 08279+5255 


3.911 


10.36 




475 


955 


394 


13.46 


8 


e 


BRI 1335-0417 


4.407 


9.77 




420 


849 


340 


13.45 


7 




BRI 0952-0115 


4.434 


9.27 




230 


436 


147 


12.38 


7 




BRI 1202-0725 


4.694 


10.16 




333 


711 


272 


13.85 


9 


f 


SDSS J1148+5251 


6.419 


9.37 




279 


579 


210 


13.43 


10 


g 



Note. — Col. (1) Quasar name. Col. (2) Redshift. We adopted the luminosity distances used by Shields et al. 2006, 
which are based on Hq = 70 km s~^ Mpc^'^, 0,^^ = 0.3, and fl\ = 0.7. Col. (3) BH mass taken from Shields et al. 2006, 
except for the following four cases (PG 0052+251, 1351+640, 1444+407, and 0157+001), whose BH masses were derived 
using Equation 2 from Shields et al., H/3 FWHM from Boroson & Green 1992, and 5100 A continuum luminosities from the 
spectrophotometry of Neugebauer et al. 1987. Col. (4) CO line width at 20% of the peak intensity. Col. (5) CO line width at 
50% of the peak intensity. Col. (6) W20 corrected for turbulence, an average inclination angle of 45°, the Tutui-Sofue effect, 
and, in the case of IF50. the luminosity dependence on line shape (see §3) using the IR luminosities from Col. (8). Col. (7) 
Stellar velocity dispersion of the bulge derived from W20 and the Vm — O"* relation. Col. (8) Total IR (8-1000 fim) luminosity. 
In the case of the high-2 quasars, they arc approximated by the far-IR luminosities taken from Solomon & Vanden Bout 2005; 
the luminosities of known lensed sources have been corrected for lensing magnification. Col. (9) Reference for the CO line 
width: (1) Schoniger & Sofue 1994; (2) Scoville et al. 2003; (3) Casoli & Loinard 2001; (4) Evans et al. 2006; (5) Planesas 
et al. 1999; (6) Weifi et al. 2003; (7) Solomon & Vanden Bout 2005; (8) Weifi et al. 2007; (9) Reichers et al. 2006; (10) 
Bertoldi et al. 2003. Col. (10) Notes: (a) I Zw 1 from Table 1; TF^'g corrected for i = 38°. (b) W20 measured manually from 
an enlarged version of the published plot; uncer tainty ~ 100 km s^^. (c) Solomon & Vanden Bout 2005 quote W50 = 280 km 
s~^ from Krips et al. 2004, but no line width information is actually given in the latter. Greve et al. 2005 attribute W50 = 
440 km s"'^ to Planesas et al. 1999 and Krips et al. 2004, but it is unclear how this value was deduced. The lensing model of 
Krips et al. suggests that the double-peaked profile of component "CO- A" arises from the disk of the host galaxy. We thus 
use the centroid velocities of the two peaks (—160 and +500 km s^^) as given in Planesas et al. 1999 to infer a full velocity 
width of 660 km s^"', which we approximate as W2n- (d) The II'so value from Weifi et al. 2003 has the smallest error bar; 
our adopted value differs slightly from the value of 400 km s^^ used by Shields et al. 2006. (e) W50 from the average value 
of 475+17 km obtained using the IRAM Plateau de Bure interferometer, which is better determined than the value from 
the IRAM 30 m. (f) VF50 differs from the value of 350 km used by Shields et al. 2006. (g) Solomon & Vanden Bout 2005 
quote W50 = 250 km s~^ from Walter et al. 2003, who fit a Gaussian profile to the combined CO(6-5) and CO(7-6) spectra 
from Bertoldi et al. 2003. The published version of Bertoldi et al.'s paper, however, actually gives W50 = 279 km s~^, which 
we adopt here. 



actual observational constraints on inclination angles, it is im- are simply biased toward more face-on orientations. To shift 
possible to disprove the prosaic alternative that these sources the ensemble of high-z points back to the local M, - cr* rela- 
tion would require the average inclination angle to be (?) = 30° 
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instead of 45°, less extreme than {i) = 15° found by Wu (2007). 



10^ 



1 1 1 1 1 1 1 1 


1 1 1 


/ : 


• 

=- 


• / 

/ 




I 


^'O / 


I 


- 

• • / 

/ 


/ / 


I 
- 



~ ^ A 

: * / A > 




- 
I 


/ / 
/ / 






r / / 
E * * / * 


• 






















^ / 
/ / A 






/ / 
/ / 






/ / 
/ / 






\ / \ / \ \ \ \ 1 


1 1 1 


1 



50 100 200 500 

a, (km s"') 



Fig. 3. — The M, — (T« relation for low-z PG quasars (triangles) and high- 
z quasars (circles). The BH masses were derived from the quasar line width 
and luminosity using the virial technique, and velocity dispersion comes from 
the CO line width. The point marked with an open circle and connected with 
a dotted line shows the new location of APM 08279+5255 with its BH mass 
adjusted to the revised luminosity lensing magnification factor of WeiB et al. 
(2007). The solid line shows the M, — cr* relation of Tremaine et al. (2002) for 
local inactive galaxies, and the dashed hiie shows an offset of 0.6 dex in mass 
from the local relation. 

Despite this impasse, the good agreement between the present 
kinematic analysis and the photometric study of Peng et al. 
(2006b) lends credence to the proposition that at high redshifts 
the growth of the BH outpaced the growth of the bulge. 

We end with a few cautionary remarks, which we hope can 
serve as a catalyst for future work. First, we should recognize 
that BH inasses derived from the virial method carry substan- 
tial uncertainty. This is especially so for masses that rely on 
the C IV A 1549 line, which might be biased with respect to 
masses obtained from H/3 (Raskin & Laor 2005). It would 
be highly desirable to acquire high-quality near-lR spectra of 
all CO-detected high-z quasars in order to reevaluate their BH 
masses using an internally consistent, homogeneous set of line 
measurements based on either H/3 or Mg II A2800 (e.g.. Earth 
et al. 2003). A number of the high-z quasars are known or 
suspected to be lensed. More accurate lensing models would 
help to bolster our confidence in their intrinsic luminosities, 
which enter crucially into the BH mass estimate via the size- 
luminosity relation. Two of the objects in Table 2 have un- 
comfortably large BH masses (M, > lO'*' Mq), larger than any 
directly measured in the local Universe. BRI 1202-0725 is 
not known to be lensed, but APM 08279+5255 is. WeiB et 
al. (2007) recently revised the luminosity magnification fac- 
tor of APM 08279+5255 from 7, the value given in Solomon 
& Vanden Bout (2005) that Shields et al. (2006) adopted, to 
100. Since the virial mass scaUng relation used by Shields et 
al. depends on L"-^, reducing the luminosity by a factor of 14 
lowers the BH mass by a factor ^ 4 (see Fig. 3). At the moment 
we do not know how many of the other lensed objects may be 
similarly affected. 

The correction for Une shape that we advocate in this study 



(Eq. 3) depends on the total IR luminosity, but we have not 
specified what actually powers this luminosity. Many IR- 
luminous galaxies are composite systems whose bolometric lu- 
minosity comes from a mixture of AGN heating and star forma- 
tion. If the empirical line shape correction is driven largely by 
the starburst component to Lir, then we may have conceivably 
overcorrected the CO line widths for the quasars, whose IR lu- 
minosity is at least partly, and some would argue predominantly 
(Ho 2005), due to AGN heating. This issue can be addressed 
by obtaining quantitative estimates of the star formation rate 
in these systems, for instance using the technique proposed by 
Ho & Keto (2007). A more straightforward remedy, of course, 
is to simply take better CO spectra so that W20 can be directly 
measured. 

Finally, the apparently excellent agreement between our re- 
sults and those of Peng et al. (2006b) is somewhat mislead- 
ing. High-redshift galaxies have steeper gravitational potential 
wells, and so a given cr* corresponds to a smaller bulge stel- 
lar mass (Robertson et al. 2006). The factor of ~4 increase in 
BH mass for high-z quasars compared to the local M, - cr* re- 
lation, therefore, actually results in an ever larger offset had the 
comparison been done in terms of stellar mass, which is the ap- 
proach taken by Peng et al. (2006b). On the other hand, Peng et 
al.'s analysis assumed a rather conservative mass-to-light ratio, 
one designed to yield the largest bulge mass. In this sense the 
factor of 4-6 increase in the BH-to-bulge mass ratio found by 
Peng et al. should be viewed as a lower Umit (C. Y. Peng 2007, 
private communications), which brings their results in closer 
agreement with ours. 



5. DISCUSSION AND SUMMARY 

The increasing realization that BHs and their energy feed- 
back play a fundamental role in the hfe-cycle of galaxies places 
great urgency in devising effective methods to diagnose high-z 
quasars and their hosts. Many approaches can be taken to tackle 
this problem, but a particularly direct one is to attempt to detect 
evolution in the well-known local scaling relations between BH 
mass and host galaxy properties. The host galaxies of high- 
z quasars, however, pose a formidable challenge for studying 
in stellar light. Not only must one contend with the enormous 
glare of the luminous AGN, but one is further handicapped by 
surface brightness dimming, band-pass shifting, and the uncer- 
tain star formation history of the galaxy. While the photometric 
properties of the stars can be studied under optimal conditions 
(e.g., in lensed systems; Peng et al. 2006b), stellar kinemati- 
cal information will be much more difficult to come by, even 
in the foreseeable future. Within this backdrop, observations of 
the large-scale cold interstellar medium of the host can serve as 
a powerful tool to probe not only its gas content but its global 
kinematics (Ho et al. 2007). From experience with galaxies 
in the nearby Universe, even a measurement as crude as an in- 
tegrated spectrum can yield an important physical parameter, 
namely the rotation velocity of the disk. The recent progress 
in the detection of CO emission in high-z quasars (Solomon 
& Vanden Bout 2005) presents an opportunity to explore this 
problem. 

Since the M, - tr* relation serves as the best local refer- 
ence benchmark, we must devise a method to relate CO line 
widths to cr* . Shields et al. (2006) tried a number of proxies, 
but were unconvinced as to which was most reUable, in large 
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part because their calibration samples were very small and be- 
cause each of the empirical correlations they employed has sig- 
nificant intrinsic scatter In the end they simply assumed cr^ 
« W5()(CO)/2.35. Wu (2007), by contrast, direcdy compared 
W5o(CO) with cr* measurements for 33 Seyfert galaxies to ar- 
rive at a conversion between the two. 

We have investigated the CO line width calibration problem 
from a very different perspective. We began by showing that 
normal galaxies obey a tight CO Tully-Fisher relation in the Ks 
band, and that we can recover the more familiar H I-based rela- 
tion provided that we use the line width measured at 20% of the 
peak (W20) instead of that measured at 50% of the peak (W50). 
The reduction in the scatter is especially impressive for IR- 
luminous galaxies, whose frequently disturbed morphologies, 
as a result of tidal interactions and mergers, produce a wide ar- 
ray of CO line shapes, and hence large differences between W20 
and W50. Indeed, we show that the line shape, as gauged by 
W20/W5Q, varies substantially and systematically with IR lumi- 
nosity. We further discuss additional corrections that need to 
be applied to the line widths in order to achieve the best pos- 
sible results. Having established W20 to be superior to W50 as 
a measure of maximum rotation velocity (v,„) of the disk, as 
well as an empirical guideline to convert W50 to W20 if the latter 
is not available, we then turn to the issue of how to relate W20 
(or equivalently v„,) to cr*. Since tr* measures the bulge and 
v„, reflects the circular velocity of the halo, it is not immedi- 
ately apparent how the two should be related. This issue was 
recently addressed by Ho (2007), who showed that although u,„ 
does not correlated as tightly with cr* as had been thought previ- 
ously, nevertheless galaxies over a wide range of Hubble types 



do follow a Um-fj* relation. This last step then completes the 
recipe to transform CO line widths to cr* . 

With these tools in hand, we revisited the CO-based M, - cr* 
relation for quasars. Low-redshift (z <^ 0.2) quasars fall roughly 
on the M, - cr* relation of inactive galaxies, albeit with consid- 
erable scatter, but high-redshift (z = 1.41-6.42) quasars statis- 
tically lie displaced from the local relation. Our finding is in 
qualitative agreement with the results of Shields et al. (2006), 
but the magnitude of the offset has been reduced with our more 
refined ct* estimates. At a given cr*, high-z quasars have BH 
masses larger by a factor of ~4 relative to local galaxies. Al- 
though we cannot rigorously rule out the possibility that we are 
being misled by a population of quasars with biased orienta- 
tions, we are encouraged by the fact that our result quantita- 
tively agrees well with that based on the photometric analysis 
presented by Peng et al. (2006b). The collective evidence sug- 
gests that early in the life-cycle of galaxies the development of 
the bulge lags behind the growth of the central BH. 
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